Analytically solvable driven time-dependent two-level quantum systems 
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Analytical solutions to the time-dependent Schrodinger equation describing a driven two-level 
system are invaluable to many areas of physics, but they are also extremely rare. Here, we present a 
simple algorithm that generates an unlimited number of exact analytical solutions. We show that a 
general single-axis driving term and its corresponding evolution operator are determined by a single 
real function which is constrained only by a certain inequality and initial conditions. Any function 
satisfying these constraints yields an exact analytical solution. We demonstrate this method by 
presenting several new exact solutions to the time-dependent Schrodinger equation. Our general 
method and many of the new solutions we present are particularly relevant to qubit control in 
quantum computing applications. 



The search for analytically solvable driven two-level 
quantum systems began shortly after the birth of quan- 
tum mechanics and continues into the present. Such sys- 
tems are ubiquitous throughout quantum physics, and it 
is notoriously difficult to acquire exact analytical solu- 
tions to the relevant time-dependent Schrodinger equa- 
tion aside from a few special cases. Perhaps the most 
famous examples of exactly soluble two-level evolution 
arc the Landau-Zener 0, Q and Rabi 0] problems. The 
latter of course has several generalizations including the 
paradigmatic Jaynes-Cummings model Q • Another well 
known exact solution discovered in the early 1930s by 
Rosen and Zener is the hyperbolic secant pulse [H[. This 
solution has proven important in the contexts of self- 
induced transparency @ and qubit control > an d nas 
since been extended to a family of analytical controls by 
a number of authors [Ioljl^ |. Some additional classes of 
soluble pulses were also discovered recent ly fiol EH ] . as 
well as a few oscillatory control examples [2l|, |22j . 

Despite these isolated successes, analytically solvable 
two-state problems have remained extremely rare. This 
fact has become more poignant in recent decades with 
the advent of quantum computation, where analytically 
solvable control pulses are especially attractive in light of 
the many advantages they offer in relation to the design 
of qubit control operations. In particular, such solutions 
can facilitate the development of controls that are both 
precise and robust without the need for long control se- 
quences [7l. lil. I231426I . However, the relatively small num- 
ber of known analytically-solvable control fields greatly 
limits the options one has when adopting an analytical 
approach to qubit gate design, and it is unlikely that any 
of the known examples will be ideal for a specific situa- 
tion, especially given that almost all of the these predate 
the inception of quantum computing. 

In this Letter, we present a completely new theoretical 
approach to the driven two-state problem. We derive an 
algorithm that produces an unbounded number of ana- 
lytically solvable two-level systems driven by a single-axis 
control field. We develop this method by showing that a 
general single-axis control field and its associated evolu- 



tion operator are both determined by a single real func- 
tion q{t), and we give the explicit functional dependencies 
on q(t). We further derive an inequality and initial condi- 
tions which q(t) must obey in order for the resulting evo- 
lution operator to be a proper solution of the Schrodinger 
equation. Any q(t) which satisfies these constraints cor- 
responds to an analytically solvable two-state problem. 
We demonstrate our method by deriving several new an- 
alytical solutions. We also determine how properties of 
q(t) translate to the control field and evolution operator. 
This 'reverse-engineering' approach is especially appro- 
priate in the context of quantum control where one typi- 
cally wishes to achieve a particular evolution by applying 
a control field whose basic features are restricted only by 
a few experimentally-imposed constraints. 

The Hamiltonian we will consider has the general form 



H 



Jit) 



(1) 



where Jit) is the control (driving) field, ft. is a constant, 
and <7 Z , o x are Pauli matrices. This Hamiltonian de- 
scribes any two-level system which is driven along a sin- 
gle axis (denoted by z) [32J. In many contexts, h can 
be interpreted as the energy splitting between the two 
levels 8, 9, 27], but in other contexts, e.g. singlet-triplet 
qubits [28143 1| . J(t) could be thought of as a time- varying 
energy splitting between the states. We parametrize the 
evolution operator corresponding to H as 

U= (« 1 ~* L )' M 9 + M a = i, (2) 

and we transform to a rotating frame in the a:-basis: 

(3) 



D ± = ±=e ±lht ' 2 iu 11 ±u 21 ). 



The functions D± then solve the following set of equa- 
tions which follow from the Schrodinger equation for the 
evolution operator U: 



Da 



(4) 



2 



These equations can be combined to yield a second-order 
differential equation for D + : 



D+ + (-ih - j/J)D+ + (J 2 /4)D + = 0. 



(5) 



At this point, one typically inserts a particular expres- 
sion for J(t) and then attempts to solve this equation 
for D + (t) to obtain the corresponding evolution opera- 
tor. Indeed, this is the manner in which most of the pre- 
viously known analytical solutions were found. In some 
cases, the precise form of J(t) was dictated by the physics 
of the problem, in other cases J{t) was chosen so that Eq. 
([5]) became a well known differential equation [33| . 

We adopt a dramatically different approach which be- 
gins by noticing that we can also view Eq. ([5]) as a differ- 
ential equation for J(t). It turns out that this equation 
can be solved exactly for arbitrary D+: 



Jit) = ± 



Da 



-iht 



\D\e~ 2lht - f /„' dt'e' 2iht ' D\(t' 



(6) 



where c is an integration constant. Given this expression 
for Jit), Eq. (j4|) then gives £)_ in terms of D + : 



2 e -2iht 



-Die 



ih 

y 



dt'e 



l„ — 2iht' 



Dun 



(7) 

Supposing the evolution begins at t = 0, we impose 
D+iO) = D-iO) = 1/V2, which in turn implies that 
c = and that we should choose the minus sign in ([7]). 
Note that imposing the initial condition at t = will not 
prevent us from obtaining solutions that span any range 
of the time domain as we will clarify later on. 

Our results so far can be interpreted as a reverse- 
engineering of the control Jit): we can choose the evolu- 
tion by picking the function D + as we like, and then use 
(|6|) and ([7J to determine -D-(i) and Jit). However, we 
must ensure that unitarity is preserved: \D + \ 2 + |Z?_ | 2 = 
1. This is automatically satisfied by the general ansatz 



D 4 



^(F-K+ht) cos$ 



D_ = e 



tK 



sin $, 



(8) 



where F, K, and $ are arbitrary real functions. These 
expressions lead to the following forms for tin and U21: 



«1 



«2 



! = ^eW-^e^cos^ + sin*), 
v2 

1 = i e i ( M /2-^( e « ? cos$ - sin$). (9) 
v2 



The initial conditions on D + , D- translate to $(0) = 
tt/4, F(0) = KiO) = 0. Eq. © imposes relations 
between F, K, and $ which can be extracted by first 
squaring both sides of this equation, differentiating the 



result with respect to time, and then equating the real 
and imaginary parts of both sides to arrive at 



F + h = if(l-tan 2 $) 
$ = iftan.Ftan$. 



(10) 



In terms of these functions, J(t) can be expressed as 

Jit) = 2XsecFtan$. (11) 

These relations further fix some additional initial condi- 
tions: $(0) = 0, F(0) = -h, and J(0) = 2Jir(0). £f(0) 
is not restricted by these relations. The next step is to 
notice that we can solve (|10[) explicitly for $ and K in 
terms of F. For $ we obtain 



sin(2$) = sec^e' l ^' it ' tanF(t ' ) , 



(12) 



where we have already chosen the integration constant 
so that the initial conditions are satisfied. One may then 
use either of the two equations in (fTU|) to solve for K. It 
then follows that once the function Fit) is specified, so 
are the control field and its evolution operator. 

As a first check we may consider the case h = 0. In this 
case, the only solution to Eq. (TT21 is to choose $ = tt/4 
and F = 0. K is then unconstrained by (TIT))) , and from 
Eq. ([9]), it is clear that we obtain a z-rotation for any 
K = i J Q dt' Jit'). We may also consider the case J = 0, 
which is realized by setting K — 0, implying that K = 0. 
It follows immediately from (TIU)) that $ = and F = —h, 
so that $ = 7r/4 and F — —ht. This solution is of course 
consistent with Eq. (| 12[) . and plugging into © reveals a 
free precession about the a;-axis as expected. 

We stress that Eq. (|12p is not simply an equation 
which gives $ once F is chosen. This equation actually 
places strong constraints on F as can be seen by noticing 
that the RHS does not generically respect the upper and 
lower bounds on sin(2$); generic choices of F will yield a 
RHS which exceeds unity. This restriction on F reflects 
the fact that our Hamiltonian, Eq. (fTJ, generates only a 
subset of all possible trajectories on the Bloch sphere. 

It is helpful to replace F by a new function q: 



(13) 



F = arctan I -p- , . 

MJ 

in which case one finds 

sin(2$) = ^q 2 + q 2 /h 2 , 



g = lhqiq + h 2 q) 



2 h 2 q 2 + q 2 



J 



q + h 2 q 



(14) 



(15) 



(16) 
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FIG. 1: (Color online) Pulse from Eq. ([20)1 with (top panel 
from left to right) a — 0, |, | and (bottom panel from left to 
right) a = 




0246 0246 



FIG. 2: (Color online) Pulse from Eq. ([22]) with (from top 
to bottom) b = — j, 0, i, 1, 2 and the corresponding evolution 
operator components. 



The initial conditions on F, K and $ translate to 

g(0) = l, 9(0) =0, q(0) = -h 2 , (17) 



and the requirement that the RHS of Eq. (|T4|) does not 
have a magnitude exceeding unity leads to 



e<h\x- q *). 



(18) 



Eqs. (fl!5]) - {P8"]) together comprise the main result of this 
paper. Any function q(t) which satisfies Eqs. (fTT]) 
and (jT8]) will produce an analytical solution to the 
Schrodinger equation, with the control field and its cor- 
responding evolution operator given by Eqs. (|13j) -(|16 |) . 
This simple prescription will enable us to generate an un- 
limited number of analytically solvable two-state prob- 
lems along with their explicit solutions. We will now 
demonstrate this by writing down several new examples. 

There is a special choice of q for which the inequality 
in Eq. (fT5)l is saturated, namely q = cos(ht); this choice 
corresponds to a pure :r-rotation with J = 0. As a first, 
more nontrivial example, we consider the choice 



q(t) = exp{-(2/a)sinh 2 (Vafti/ 2 )} , 



(19) 



where a is any real number satisfying a < 2. It can be 
verified that {19]) satisfies Eqs. {T7]) and {18]). Eq. {IS]) 
then gives 



J(t) 



h [(1/a) sinh 2 ( v / a^) - 2 sinh 2 ( v / aW/ 2 )] 

y / e (4/a)sin h 2( v ^/ lt /2) _ s i n h 2 (^/Eht) 



(20) 

and this result is plotted for different values of a in Fig. [T] 
The corresponding evolution operator can be computed 
straightforwardly from Eqs. dT3^) - (IT51) . An interesting fea- 
ture of this solution is that it describes a single pulse 



when a > and an oscillatory control field when a < 0. 
This solution illustrates that the complexity of J(t) and 
U(t) tends to be comparable to that of q(t), in contrast 
to most of the analytical solutions known prior to this 
work, where simple pulse shapes often yield evolutions 
governed by special functions. 

A second new nontrivial example arises from the choice 



«(*) = 



1 



1 + 6 L 



bcos(ht)], 



which upon inserting into Eq. (fT6l) yields 



J(t) 



h 3 t 2 e -hV/2 



y/l-(l + hH 2 )e- hH * + 2b X {t) ' 



(21) 



(22) 



with x(t) 



-h"t 



!/2 [cos(M) + htsm(ht)]. Plots of 



this pulse and its corresponding evolution operator for 
different values of b are shown in Fig. [5] where it is ap- 
parent that these are smooth, Gaussian-like pulses with 
b controlling the magnitude and width. These pulses of- 
fer an attractive alternative to the widely used hyper- 
bolic secant pulses [B| : they converge to zero much more 
rapidly and are thus more localized, and the correspond- 
ing evolution operator is expressed in terms of elemen- 
tary functions rather than the hypergeometric function. 
Gaussian pulses have proven useful in the context of error 
suppression in superconducting qubit gate design 24 - 2(| , 
however these were not analytically solvable and required 
numerical optimization techniques. 

The two examples we have given so far already illus- 
trate that properties of q(t) are reflected in J(t). In the 
context of quantum control, experimental constraints of- 
ten require J(t) to be smooth, bounded, and well local- 
ized in time. Smoothness of J(t) is guaranteed by choos- 
ing a smooth q(t). In addition, it is clear from Eq. (|16[) 



that a given q(t) will produce a well defined pulse (i.e. 
J(t) — ¥ as t — > ±00) if q and g — >• as i — ► ±00, 
or if q — > —h 2 q in this limit (the former condition is 
satisfied by Eq. $T3§ and the latter by Eq. (|2"T])l The 
strict inequality < q 2 < h 2 (l — q 2 ) for t > further 
ensures that J(£) will be a bounded function. In some 
cases (e.g. singlet-triplet qubits [3l|), one also needs to 
impose positive control, J(t) > 0; as seen from Eq. (fl"6|) . 
this requires q > -h 2 q. Eq. (19]) with a > and Eq. (|2Tj) 
satisfy the above criteria, in which case we have a single, 
positive, bounded pulse as confirmed by the top panels 
of Fig. [1] and the top left panel of Fig. [2] When a < 
in Eq. (fl9|) . q(t ) = for some to > 0, leading to the 
periodic behavior shown in the lower panels of Fig. [TJ 

The examples of Eqs. (|19[) and ([21]) also exhibit the 
connection between q(t) and the corresponding evolu- 
tion operator. For instance, these examples contain 
even control functions, J(—t) = J(t), which follows di- 
rectly from the evenness of the chosen q(t). For such 
pulses, the evolution operator given by Eqs. ((2|) and (|9]) 
describes half the evolution due to the pulse since we 
have imposed f (0) = 1. If we wish to instead impose 
U(—tf) — 1 for some tf > 0, then the full evolution op- 
erator describing the evolution from t = —tf onward is 
U tot (t) = U(t)W(-t f ). At t = t f , one then finds that 
T r{U tot {tf)a y } = a nd Tr{U tot (t f )a z } ~ sra(2$(t/)) = 
y/q 2 (tf) + q 2 (tf)/h 2 , implying that the pulse effects a 
rotation about an axis in the x — z plane which de- 
pends on the behavior of q at t = tf. For well de- 
fined pulses in particular, if tf is sufficiently large and 
q — > Acos(ht) + Bsin(ht) as t — > 00, then the rotation 
axis is determined by A and B, where A = B = yields 
an x-rotation. It is then clear from Eq. (fl9| that the 
pulse family of Eq. (|20[) implements x-rotations for any 
a > 0, while Eq. (|21l) reveals that the pulses given in Eq. 
(122|l implement rotations about various axes in the x — z 
plane depending on b. In particular, denoting the rota- 
tion axis by (n x , 0, n z ) and the angle by 9, the lower left 
panel of Fig. [2]reveals that lm[U to t, 11 (tf)] = —n z sm(9/2) 
quickly saturates to a constant ^-dependent value beyond 
htf > 3. Given n z and 0, one can first choose b to fix 
the combination n z sm(9/2) and then tune tf > 3//i to 
achieve the target rotation. 

We can systematically find q's which satisfy Eq. (TT5)) 
by first choosing a function V(q) such that < P(q) < 
1 — q 2 and then solving q 2 = h 2, P(q). Since this equation 
is homogeneous in q, it can be integrated directly: 



lit 



1 dq' 



= W(q). 



(23) 



The fact that V(q) is strictly nonnegative guarantees that 



nl dq' 

W(q) > I — ; = arccosg. 



.. , (24) 

A 7 ? 

The initial conditions on q(t) become the condition 




FIG. 3: (Color online) Pulse from Eq. |g7]| with (left panel, 
from top to bottom) a = 2v2, 2, \/2, 1, l/y/2 and (right panel, 
from top to bottom) a = l/y/2, 0.6, 0.5, 0.4, 0.3. 
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FIG. 4: (Color online) Control from Eq. J28]) with (left) a 
0.1 and (right) a = 0.5. 



W(q) — > \J2 — 2q as q — > 1. Any invertible function 
W(q) satisfying this boundary condition and the inequal- 
ity ([24]) automatically produces an analytical solution to 
the Schrodinger equation with q(t) = W^ 1 (ht). 

To give an example using this approach, we choose 



W(q) = (l/a)arctanh(aA/2 - 2q), (25) 
where a is a real constant. Upon inverting, we obtain 



q(t) =1- — t&nh 2 (aht), 
2a l 



(26) 



h [14a 2 



(2a 2 - l)cosh(2aW)] sech 2 (aht) 



J(t) = 

2coth(aht)Jia 2 [l - sech 4 (aht)] - tanh 2 (aht) 

(27) 

The corresponding control pulses are shown in Fig. [3] 
The pulse with a = l/y/2 asymptotes to zero at large 
times, while pulses with a > l/y/2 asymptote to the pos- 
itive constant (2a 2 — 1)/V4a 2 — 1. These latter pulses 
could be relevant in experimental situations in which it 
is not possible to turn off J(t) completely, as in the case 
of singlet-triplet qubits 28|-3(|- Pulses with a < l/y/2 
are defined on a strictly finite time interval as shown in 
the right panel of Fig. [3] 

As a final example, we make the choice 

W(q) = cos -1 {l - (a+ 1/a) [tan _1 (ag) - tan -1 a] } , 

yielding 

q(t) = (1/a) tan {tan -1 a - (2a/(l + a 2 )) sin 2 (ht/2)} . 

(28) 
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The resulting J(t) is a straightforward combination of el- 
ementary functions that is plotted in Fig. |4] where it is 
apparent that this control resembles two superimposed 
sinusoids. Varying the parameter a changes the ampli- 
tude and average value of these sinusoids, but not their 
frequency. This type of solution is relevant, e.g., for de- 
termining the properties of two-level fluctuators which 
cause decoherence in superconducting qubits [HI, [53] . 

In conclusion, we have shown how to systematically 
obtain an unlimited number of analytically solvable con- 
trols and have provided explicit analytical formulas for 
their corresponding evolution operators. This vastly in- 
creases the number of known analytical solutions to the 
two-state Schrodinger equation, providing a powerful tool 
in robust quantum gate design and in the myriad of other 
physical contexts in which the two-state problem arises. 

This work is supported by LPS-NSA and IARPA. 
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